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qj ' Abstract 

Evolution of the concept known in the theoretical physics as the Renormalization 
Group (RG) is presented. The corresponding symmetry, that has been first intro- 
duced in QFT in mid-fifties, is a continuous symmetry of a solution with respect to 
transformation involving parameters (e.g., of boundary condition) specifying some 
particular solution. 

After short detour into Wilson's discrete semi-group, we follow the expansion 
of QFT RG and argue that the underlying transformation, being considered as a 
reparameterisation one, is closely related to the self-similarity property. It can be 
treated as its generalization, the Functional Self-similarity (FS). 

Then, we review the essential progress during the last decade of the FS concept 
in application to boundary value problem formulated in terms of differential equations. 
A summary of a regular approach recently devised for discovering the RG = FS sym- 
metries with the help of the modern Lie group analysis and some of its applications 
are given. 

As a main physical illustration, we give application of new approach to solution 
for a problem of self-focusing laser beam in a non-linear medium. 
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1 The Bogoliubov renormalization group 
1.1 Historical introduction 

1.1.1 The discovery of the renormalization group 

In 1952-1953 Stiickelberg and Petermann [jlj discovered^ a group of infinitesimal trans- 
formations related to finite arbitrariness arising in S'-matrix elements upon elimination 
of ultraviolet (UV) divergences. These authors used the notion of normalization group as 
a Lie transformation group generated by differential operators connected with renormal- 
ization of the coupling constant e. 

In the following year, on the basis of (infinite) Dyson's renormalization transformations 
formulated in the regularized form, Gell-Mann and Low 0] derived functional equations 
(FEs) for the QED propagators in the UV limit. The appendix to this article contained 
the general solution (obtained by T.D. Lee) of FE for the renormalized transverse photon 
propagator amplitude d(Q 2 /X 2 ,e 2 ) (A - cutoff defined as a normalization momentum). 
This solution was used for a qualitative analysis of the quantum electromagnetic inter- 
action behaviour at small distances. Two possibilities, namely, infinite and finite charge 
renormalizations were pointed out. However, paper || paid no attention to the group char- 
acter of the analysis and of the qualitative results obtained there. The authors missed a 
chance to establish a connection between their results and the QED perturbation theory 
and did not discuss the possibility that a ghost pole solution might exist. 

The decisive step was made by Bogoliubov and the present author [|], [|, [J in 1955 []. 
Using the group properties of finite Dyson transformations for the coupling constant, 
fields and Green functions, they derived functional group equations for the renormalized 
propagators and vertices in QED in the general (i.e., with the electron mass taken into 
account) case. 

In the modern notation, the first equation 

_(x y \ Q 2 m 2 

-,-;«(t,,;a)j ; * = -, y=- 



a(x, y; a) = a(j,~; a(t, y;a)) ; x = — , y = — , (1) 



is that for the invariant charge (now widely known also as an effective or running coupling) 
a = ad(x, y; a = e 2 ) and the second — 

s(x, y\ a) = s(t, y; a) s (-, -; a(t, y; a)) (2) 

— for the electron propagator amplitude. 

These equations obey a remarkable property: the product e 2 d = a of the electron 
charge squared and the photon transverse propagator amplitude enters into both FEs. 



1 For a more detailed exposition of the RG early history, see our recent reviews 
2 See also two survey papers published in English in 1956. 



This product is invariant with respect to finite Dyson's transformation (as it is stated by 
Eq. ([!])) which now can be written in the form 

Rt : {/i 2 -> t/i 2 ,a — ► a(t, y;a)j . (3) 

We called this product invariant charge and introduced the term renormalization 
group. 

Let us emphasize that, unlike in Refs.[|l[ ||, in the Bogoliubov formulation there is no 
reference to UV divergences and their subtraction or regularization. At the same time, 
technically, there is no simplification due to the massless nature of the UV asymptotics. 
Here, the homogeneity of the transfer momentum scale Q is explicitly violated by the 
mass m. Nevertheless, the symmetry with respect to transformation R t (even though a 
bit more involved) underlying RG is formulated as an exact property of the solution. This 
is what we mean when using the term Bogoliubov renormalization group or renormgroup 
for short. 

The differential Lie equations for a and for the electron propagator 

da{x,y;a) (y \ ds(x,y;a) (y \ 

— — =p [-, a{x,y a)) — — = 7 -, a{x,y a) )s(x,y a) , (4) 

omx \x J omx \x J 

with 

/3{y,a) = — , 7U/, a) = ^7 at £ = 1 (5) 

were first derived in [|J by differentiating FEs ([TJ) and (0) over x at the point t = x. On 
the other hand, by differentiating the same equations over t one obtains || 

Xa(x,y;a) =0 ; X s(x,y; a) = 7(2/, a)s(x,y; a) (6) 

with 

X = xd x + yd y - (3(y, a)d a (d x = d/dx) , (7) 

the Lie infinitesimal operator. 

1.1.2 Creation of the RG method 

Another important achievement of || consisted in formulating a simple algorithm for 
improving an approximate perturbative solution by combining it with Lie group equations 
- for detail, see below Section |1.3j 

In our adjacent publication this algorithm was effectively used to analyse the UV 
and infrared (IR) behaviour in QED. In particular, the one- loop UV asymptotics of the 
photon propagator as well as the IR behavior of the electron propagator in the transverse 
gauge 

a g)(x;a) = T — |^- , s(x,y;a)^(p 2 /m 2 -iy^ (8) 

3tt 



were derived. At that time, these expressions, summing the leading log's terms were al- 
ready known from papers by Landau with collaborators [|§. However, Landau's approach 
did not provide a means for constructing subsequent approximations. 

The simple technique for calculating higher approximations was found only within the 
new renormgroup method. In the same paper, starting with the next order perturbation 
expression a^ t (x; a) containing the a 3 lnx term, we arrived at the second renormgroup 
approximation (see below Section |1.3.2|) 



a 



" rg(a:;a) 1-f lnx + |fln(l-f lnx) (9) 

which performs infinite summation of the a 2 (aln) n terms. This two-loop solution for 
invariant coupling first obtained in || contains the nontrivial log-of-log dependence which 
is now widely known of the "next-to-leading logs" approximation for the running coupling 
in quantum chromodynamics (QCD) — see, below, Eq. fl2T|) . 

Comparing (|9[) with (§), one concludes that two- loop correction is essential in the 
vicinity of the ghost pole at x\ = exp (3ir/a). This also shows that the RG method is 
a regular procedure, within which it is easy to estimate the range of applicability of its 
results. 

Quite soon, this approach was formulated || for the case of QFT with two coupling 
constants. To the system of FEs for two invariant couplings there corresponds a coupled 
system of nonlinear differential equations (DEs). The last was used |]10[ to study the UV 
behavior of the n — N interaction at the one-loop level. 

Thus, in Refs. ||, |], ||] an d ]1(| RG was directly connected with practical computations 



of the UV and IR asymptotics. Since then, this technique, the renormalization group 
method (RGM)Q, has become the sole means of asymptotic analysis in local QFT. 

1.2 The Bogoliubov RG: Symmetry of a solution 

The RG transformation. Generally, RG can be defined as a continuous one-parameter 
group of specific transformations of a partial solution (or the solution characteristic) of 
a problem, a solution that is fixed by a boundary condition. The RG transformation 
involves boundary condition parameters and corresponds to some change in the way of 
imposing this condition. 

For illustration, imagine one- argument solution characteristic f(x) that has to be 
specified by the boundary condition f(x ) = f . Formally, represent a given characteristic 
of a partial solution as a function of boundary parameters as well: f(x) = f(x,x ,f ). 
This step can be treated as an embedding operation. Without loss of generality / can be 
written in a form of a two-argument function F(x/xo, fo) with the property -F(l, 7) = 7. 

The RG transformation then corresponds to a changeover of the way of parameteriza- 
tion, say from {x , fo} to {xi, /1} for the same solution. In other words, the x argument 

3 Being summarized in the special chapter of the first edition of monograph JEll. 



value, at which the boundary condition is given, can be changed for X\ with f(x±) = f\. 
The equality F(x/x , f ) = F(x/xx, /i) now reflects the fact that under such a change 
the form of the function F itself is not modified. Noting that f\ = F(xi/xq, /q) , we get 

F(t, fo) = F(£/t, F(t, /„)) ; e = x/x , t = Xl /x . 

The group transformation here is { £ — > £/£, /o — > -F(t, /o) } • 

The renormgroup transformation for a given solution of some physical problem in the 
simplest case can now be defined as 

a simultaneous one-parameter transformation of two variables, say x and g, by 

R t : { x -»• x' = x/t , -> flr' = #(£, #) } , (10) 

the first being a scaling of a coordinate x (or reference point) and the second — a more 
complicated functional transformation of the solution characteristic. The equation 

g(x,g) =g(x/t, g(t,g)) (11) 

for the transformation function ~g provides the group property T Tt = T T T t of the trans- 
formation flip]) . 

They are just the RG FEs and transformation for a massless QFT model with one 
coupling constant g. In that case x = Q 2 / fj, 2 is the ratio of a 4-momentum Q squared to 
a "normalization" momentum jjl squared and g, the coupling constant. 

The RG transformation (|T0|) of a QFT amplitude s is of the form (compare with 
Eq.©) 

Rt -s(x,g) = e- lntx s(x,g) = s(x/t, g(t,g)) = z~ l s(x,g); z s = s(t,g). (12) 

Several generalizations are in order. 

a. "Massive" case. For example, in QFT, if we do not neglect mass m of a particle, we 
have to insert an additional dimensionless argument into the invariant coupling g which 
now has to be considered as a function of three variables: x = Q 2 /fi 2 , y = m 2 /fi 2 , and 
g. The presence of a new "mass" argument y modifies the group transformation ( |10D and 
the FE flnD 

y = -, g' = g{t,y;g) j ; g{x,y;g) = g(j,-;g{t,y;g)). (13) 



/?,: {x' = - t 



Here, it is important that the new parameter y (which, physically, should be close to the 
x variable, as it scales similarly) enters also into the transformation law of g . 

If the considered QFT model, like QCD, contains several masses, there will be several 
mass arguments y ->• {y} = y t , y 2 , . . . , y n ■ 



b. Multi-coupling case. A more involved generalization corresponds to transition to the 
case with several coupling constants: g —>■ {g} — gx, . . .,gj. ■ Here, there arise a "family" 
of effective couplings 



i — 1, 2, . . . , k , 



(14) 



9 ^ {g} , 9i = 9i(x,y;{g}) ; 
satisfying the system of coupled functional equations 

9i{x, y; {g}) = g t (x/t , y/t ; { g(t, y; {g}) } ) . (15) 

The RG transformation now is 

Rt : {x->x/t, y^y/t, {g} -► {g(t)} } , 9i (t) = gi(t,y;{g}) . (16) 



1.3 Renorm-group method 

1.3.1 The algorithm 

The idea of the approximate solution marriage j|, [5| to group symmetry can be realized 
with the help of group DEs. If we define f3 and 7 (the so-called "generators" on physical 
slang) from some approximate solutions and then solve evolutional DEs, we obtain the 
RG improved solutions that obey the group symmetry and correspond to the approximate 
solutions used as an input. 

Now we can formulate an algorithm of improving an approximate solution. The proce- 
dure is given by the following prescription which we illustrate by a massless one-coupling 
case (|) and (g): 

Assume some approximate solution g appr (x,g) ,s appr (x,g) is known. 

1. On the basis of Eq.([5|) define the beta- and gamma-functions 



P(g) = ^70appr(f.0) 



€=1 



I \ dof ^ It 

i\g) = ^ appr (£, 



5=1 



(17) 



2. Integrate the first of Eqs.(f|), i.e., construct the function 

dcf f 9 dr 



fig) 



P{r) 



3. Resolve the obtained equation, i.e., 

g rK (x,g) = f' 1 {f(g) + ^x} . 



(19) 



4. Integrate the second of Eqs.(§) using this expression g rg in its r.h.s. to obtain 
s rg (x,g) in the explicit form. 

5. Then, the expressions g rg and s rg precisely satisfy the RG symmetry, i.e., they are 
exact solutions of Eqs. ([[]]) and ([12]) corresponding to (? appr and s appr used as an input. 



1.3.2 Simple illustration 



ttW 



appr 



For illustration, take the simplest perturbative expression g pt — g — g f3\ In x for g 

and Sp t = 1 — g7ilnx. Here, (3(g) = —fag 2 , 7(g) = —715 and integration of (||) gives 
explicit expressions 

9 { rg(x,g) = — — , 8®(x,g) = (g(x,g)/g) Vl ] v 1 =j 1 //3 1 , (20) 

which, on the one hand, exactly satisfy the RG symmetry and, on the other, being ex- 
panded in powers of g, correlate with g pt and s pt . 

Now, on the basis of geometric progression (f20|), let us present the 2- loop perturbative 
approximation for g in the form g pt = g — g 2 j3i In x + g 3 ((3 2 In x — fa In x) . By using this 
expression as an input in Eq.([L7|), we have P^(g) = —Pig 2 — fag 3 and then (step 2) 

To make the last step, we have to start with the equation 

/(%g)(x,s)] = /( 2 )( 5 ) + /Mn* 

which is a transcendental one and has no simple explicit solution^]. Due to this, one usually 
resolves this relation approximately. Take into account that the second, logarithmic, 
contribution to p 2 \z) is a small correction to the first one at bz <C 1. Under this 
reservation, we can substitute the one-loop RG expression (p0|) instead of g^ 2 ) into this 
correction and obtain the explicit "iterative" solution 

77(2) = i. ■ 7 = lnr (211 

*k 1 + 0ftl + g(fa/fa) In [1 + gfal] ' • l J 

An analogous procedure for s pt = 1 — (771 In a; + g 2 (71(71 + l)/3f In 2 x — 72 In 2 J yields 



( ,, g^g|W)) .., Q , s ^ „ 29 , A72- 7 i& 
g = 5f ) — ^ = 9 " 2 = tf 2 " 



S W = v g-l- y with %) = ^e" 29 and z/ 2 = ^ ,z ^ '^ . (22) 

These results are interesting from several aspects. 



First, being expanded in g and gl powers, they produce an infinite series containing 
"leading" and "next-to-leading" UV logarithmic contributions. 

Second, they contain anontrivial analytic dependence ln(l+g/3il) ~ ln(ln<5 2 ) which 
is absent in the perturbation input. 

Third, being compared with the one- loop solution, Eq. (|20|) , they demonstrate an 
algorithm of subsequent improving of accuracy, i.e., of RGM regularity. 



4 It can be expressed via special, Lambert, VF-function: W(z) exp w ^ = z; see, e.g., Rcf. [p~2 



1.3.3 RGM usage in QFT 

As we have seen, QFT perturbation expression of finite order does not obey the RG sym- 
metry. On the other hand, it was shown that the one-loop and two-loop approximations, 
used as an input for the construction of "generators" /3(g) and j(g), yield expressions 
( p0|) , fl2lD and (|2|) that obeys the group symmetry and exactly satisfy FEs (|ll|) and ([12]) . 
More generally, one can state the following logical structure of the RGM procedure. 

- Solving group equation(s) for invariant coupling(s) g rg (x, g) using some approximate 
solution g pt as an input. 

- Obtaining RG solutions for some other QFT objects (like vertices and propagator 
amplitudes) on the basis of the expression(s) for g Tg just derived. Typically, they satisfy 
the equation 

XM(x,y,g)=j(y,g)M(x,y,g). (23) 

General structure of the corresponding solutions has the form 

M(x,y;g) = z£(y,g) M(x/y , g(x,y;g)). (24) 

Note that the function M. in the r.h.s. depends only on the RG invariants, that is 
on the first integrals of the RG operator X introduced in Eqs.(^J) and (0). It satisfies 
homogeneous partial differential equations (PDEs) X M. = 0. For the RG invariant 
objects, like observables, Zm — 1, 7 = 0. 

Now we can resume the RGM properties. The RGM is a regular procedure of combining 
dynamical information (taken from an approximate solution) with the RG symmetry. The 
essence of RGM is the following: 

1) The mathematical tool used in RGM is Lie differential equations. 

2) The key element of RGM is possibility of an (approximate) determination of "genera- 
tors", like (3(g), j(g), from dynamics. 

3) The RGM works effectively in the case when the solution has a singular behaviour. It 
restores the structure of singularity compatible with the RG symmetry. 

2 Evolution of Renormalization Group 

In the 70s and 80s RG ideas have been applied to critical phenomena: spontaneous mag- 
netization, polymerization, percolation, non-coherent radiation transfer, dynamic chaos, 
and so on. Less sophisticated motivation by Wilson in spin lattice phenomena (than in 
QFT) made this "explosion" of RG applications possible. 

2.1 Renormalization Group Evolving 



2.1.1 Kadanoff- Wilson RG in critical phenomena 

a. Spin lattice. The so-called renormalization group in critical phenomena is based on the 



Kadanoff- Wilson procedure [13], |14| referred to as "decimation" or "blocking". Initially, 
it emerged from the problem of spin lattice. 

Imagine a regular (two- or three-dimensional) lattice consisting of N d , d = 2, 3 cites 
with an 'elementary step' a between them. Suppose that at every site a spin vector a is 
located. The Hamiltonian, describing the spin interaction between nearest neighbours 

H = kj^a-i- <Xj±i 

i 

contains k, the coupling constant. A statistical sum is obtained from the partition func- 
tion, S =< exp(—H/9) > aver • 

To realize blocking, one has to perform the "spin averaging" over block consisting of 
n d elementary sites. This step diminishes the number of degree of freedom from N d to 
(N/n) d . It also destroys the small-range properties of a system, in the averaging course 
some information being lost. However, the long-range physics (like correlation length 
essential for phase transition) is not affected by it, and we gain simplification of the 
problem. 

As a result of this blocking procedure, new effective spins S arise in new sites forming 
a new effective lattice with a step na. We arrive also at the new effective Hamiltonian 

H e s — K n }] Sj • S/ ± i + AH , 

with the effective coupling K n between new spins S/ of new neighbouring sites; K n has to 
be defined by the averaging process as a function of k and n . Here, AH contains quartic 
and higher spin forms which are irrelevant for the IR (long-distance) properties. Due to 
this, one can drop AH and conclude that the spin averaging leads to an approximate 
transformation, 

i I 

or, taking into account the "elementary step" change, to {a — > na, k — > K n }. The latter 
is the Kadanoff- Wilson transformation. It is convenient to write down the new coupling 
K n in the form K n = K(l/n, K). Then, the KW transformation reads 

KW n : {a^na, k -> K n = K (1/n, k)} . (25) 

These transformations obey composition law KW n ■ KW m = KW nm if the relation 

K(x,k) = K(x/t,K(t,k)), x = l/nm, t = 1/n . (26) 

holds. This is very close to RG symmetry. 
We observe the following points: 

10 



• The RG symmetry is approximate (due to neglecting AH). 

• The transformations KW n are discrete. 

• There exist no reverse transformation to KW n . 

• Transformations KW n relate different auxiliary models. 

Hence, the 'Kadanoff-Wilson renormalization group' (KW-RG) is an approximate and 
discrete semi-group. For a long-distance (IR limit) physics, however, A(l/n) is small and 
it is possible to use differential Lie equations^]. 

b. Polymer theory. In polymer physics, one considers statistical properties of polymer 
macromolecules which can be imagined as a very long chain of identical elements (with 
the number of elements N as big as 10 5 ). Molecules are swimming in a solvent and form 
globular s. This big molecular chain forms a specific pattern resembling the pattern of a 
random walk. The central problem of the polymer theory is very close to that of a random 
walk and can be formulated as follows. 

For a long chain of iV "steps" (the size of step = a), one has to find the "chain size" 
Rn, the distance between the "start" and the "finish" points (the size a of globular), with 
the distribution function /(</») of angles between the neighboring elements being given. 

For large N values, the molecular size Rn follows the power Fleury law Rn ~ N u 
with u, the Fleury index. When N is given, Rn is a functional of /(</>) which depends 
on external conditions (e.g., temperature T, properties of solvent, etc. ). If T grows, Rn 
increases and at some moment globulars touch one another. This is the polymerization 
process which is very similar to a phase transition phenomenon. 

The Kadanoff-Wilson blocking ideology has been introduced in physics of polymers 



by De Gennes [15|. The key idea is a grouping of n neigbouring elements of a chain into 
a new "elementary block". It leads to the transformation {1 — > n ; a — > A n } which is 
analogous to one for the spin lattice decimation. This transformation must be specified 
by a direct calculation which gives an explicit form of A n = a(n,a). Here, we have a 
discrete semi-group. Then, by using the KW-RG technique, one finds the fixed point, 
obtains the Fleury power law and can calculate its index v. 

The essential feature of a polymer chain is the impossibility of a self-intersection. 
This is known as an excluded volume effect in the random walk problem. Generally, the 
excluded volume effect yields some complications. However, inside the alternative, the 



QFT RG approach to polymers |16[, it can be treated rather simply by introducing an 



one more argument which is similar to finite length L in the transfer problem or particle 
mass m in QFT. 

Besides polymers, the KW-RG technique has been used in some fields of physics, like 
percolation, non-coherent radiation transfer |l7j], dynamical chaos [18| and some others. 



5 In application of these transformations to critical phenomena, the notion of a fixed point is important. 
Generally, a fixed point is associated with power-type asymptotic behavior. Note here that, contrary to 



the QFT case considered in Section 1.3.2, in phase transitions we deal with the IR stable point. 

11 



2.1.2 Bogoliubov symmetry outside QFT 

Meanwhile, the original QFT-RG approach proliferated into some other parts of theoret- 
ical physics. In the late 50s, it was used |19| for summation of Coulomb singularities in 



Bogoliubov's theory of superconductivity based on the Frohlich electron-phonon interac- 
tion. Twenty years later it was used in the theory of turbulence. 

a. Turbulence. To formulate the turbulence problem in terms of RG, one has to perform 
the following steps p0| , |2Tf : 



1. Introduce the generating functional for correlation functions. 

2. Write down the path integral representation for this functional. 

3. By changing the functional integration variable, find the equivalence of the statistical 
system to some quantum field theory model. 

4. Construct the system of Schwinger-Dyson equations for this equivalent QFT model. 

5. Perform the finite renormalization procedure and derive the RG equations. 

Here, the reparameterization degree of freedom, physically corresponds to a change of 
long wave-length cutoff which is built-in into the definition of a few effective parameters. 

b. Weak shock wave. Another example can be taken from hydrodynamics. Consider 
a weak shock wave in the one-dimensional case of a large distance I from the starting 
(implosion) point. The dependence of velocity v of a matter as a function of I at a given 
moment of time t has a simple triangular shape and can be described by the expression 

v (l) = -V at I < L ; = for I > L , 

where L = L(t) is the front position and V = v(L) - the front velocity. They are functions 
of time. In the absence of viscosity, the "conservation law" LV = Const, holds. Due 
to this, they can be treated as functions of the front wave position L = x, V = V(x) 
as well. If the physical situation is homogeneous, then the front velocity V(x) should 
be considered as a function of only two additional relevant arguments - its own value 
Vo = V(x ) at some precedent point (x < x) and of the x coordinate. In can be written 
down in the form : V(x) = G(x/x , Vo) • If we pick up three points x , x\ and x 2 (for 
details, see Refs.[^, p3|), then the initial condition may be given either at x or x±. Thus, 
we obtain the FE equation equivalent to (|TT1) 

V 2 = G(x 2 /x , V ) = G(x 2 /x 1 , Vi) = G(x 2 /x 1 , G(x 1 /x , V )) . 



c. One-dimensional transfer. A similar argument has been done by Mnatzakanian [26 



in the transfer problem at one dimension. Imagine a half-space filled with a homogeneous 

12 



medium on the surface of which some flow (of radiation or particles) with intensity go falls 
from the vacuum half-space. 

Follow the flow as it moves inwards the medium at the distance I from the boundary. 
Due to homogeneity along the I coordinate, the intensity of the penetrated flow g{l) 
depends on two essential arguments, g{l) = G(l,g ). The values of the flow at three 
different points go (on the boundary), g\ and #2 > 9\ can be connected with each other 
by the transitivity relations, g\ = G(X,go) , gi = G(X + I, go) = G(l,gi) , which lead to 
the FE 

G(l,g) = G(l-X,G(X,g)). (27) 

Performing a logarithmic change of variables I = lnx, A = hit, G(l,g) = g(x,g) , we 
see that ( PT| ) is equivalent to (JTTj) . 

Consider now intensity of a reverse flow, that is total amount of particles at the point 
I moving in the backward direction. It is completely defined by go an d can be written 
down as R(l,g ). This function can be represented in the form R(l,g) = Ro(g)N(l,g) 
with R = R(0,g) and function N "normalized" on the boundary N(0,g) = 1. Playing 
the same game with transitivity, we arrive at FE 

N(l, g) = Z(l, g)N(l - A, G(l, G(X, g)) ; Z = Ro{gi)/Ro(g) (28) 



related to Eq.(|T2D by logarithmic change of variables. One can refer to (^7]) and (^g) as 
to the additive version of RG FEs and to previous equations of Section 1, like (|TT|), flT2| ) 
and (p~3| ) as to the multiplicative one. 

The transfer problem admits a modification connected with discrete inhomogeneity: 
imagine the case of two different kinds of homogeneous materials separated by the inner 
boundary surface at I = L. The point of breaking / = L may correspond to the boundary 
with empty space, and resulting equation is equivalent to Eq. ([l3|) . 

One more generalization is related to "multiplication" of argument g as expressed by 
Eq . (|14D . Physically, this relates to the case of radiation on different frequencies u>i , i = 
1, 2, . . . k (or particles of different energies or of different types). 

Take the case of k = 2 and suppose that the material of the medium has such prop- 
erties that the transfer processes of the two flows are not independent. In this case, the 
characteristic functions of these flows G and H are dependent on both the boundary 
values go and h and can be taken as functions g{l) = G(l,go, h ) , h{l) = H(l,g , ho) . 

After a group operation I —> I — A , we arrive at a coupled set of functional equations 

G(/+A, g, h) = G{1, g x , h x ) , H(l+X, g, h) = H(l, g x , h x ) ; g x = G{X, g, h) , h x = H(X, g, h) 

which is just an additive version of system (|15D at k = 2. 

Now we can make the important conclusion that a common property yielding func- 
tional group equations is just the transitivity property of some physical quantity with 
respect to the way of giving its boundary or initial value. 

Hence, the RG symmetry is not a symmetry of equations but a symmetry of solution, 
that is of equations and boundary conditions considered as a whole. 
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2.2 Difference between Bogoliubov RG and KW— RG 

As we have mentioned above, the RG ideas expanded in diverse fields of physics in two 
different ways: 

• via direct analogy with the Kadanoff-Wilson construction (averaging over some 
set of degrees of freedom) in polymers, non-coherent transfer and percolation, i.e., 
constructing a set of models for a given physical problem. 

• via finding an exact RG symmetry by proof of the equivalence with a QFT model 
(e.g., in turbulence |2(| |21|), plasma turbulence [p7| ) or by some other reasoning 



(like in a transfer problem). 

To the question Are there different renormalization groups? the answer is positive: 

1. In QFT and some simple macroscopic examples, RG symmetry is an exact sym- 
metry of the solution formulated in its natural variables. 

2. In turbulence, continuous spin-field models and some others, it is a symmetry of an 
equivalent QFT model. 

3. In polymers, percolation, etc. , (with KW blocking), the RG transformation is 
a transformation between different auxiliary models (specially constructed for this 
purpose) of a given system. 

As we have shown, there is no essential difference in the mathematical formulation. 
There exists, however, a profound difference in physics: 

- In the cases 1 and 2 (as well as in some macroscopic examples), the RG is an exact 
symmetry of a solution. 

- In the Kadanoff-Wilson type problem (spin lattice, polymers, etc. ), one has to 
construct a set M. of models Mj. The KW-RG transformation 

KW n Mi = M ni , with integer n (29) 

is acting inside a set of models. 

2.3 Functional self-similarity 

The RG transformations have close connection with the concept of self-similarity The 
self-similarity transformations for problems formulated by nonlinear PDEs are well known 
since the last century, mainly in dynamics of liquids and gases. They are one parameter 
A transformations defined as a simultaneous power scaling of independent variables z = 
{x, t, . . .} , solutions fk(z) and other functions Vi(z) (like external force) 

S x : {x' = x\,t' = t\\ f' k = X^f k , V! = \ Vi Vi } 
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entering into the equations. 

To emphasize their power structure, we use a term power self-similarity = PS. Accord- 



ing to Zel'dovich and Barenblatt, |25|] the PS can be classified as: 

a/ PS of the 1st kind with all indices a, ...cp, u, ... being integers or rational (Rational 
PS) that are usually found from the theory of dimensions; 

b/ PS of the 2nd kind with irrational indices (Fractal PS) which should be defined 
from dynamics. 

To relate RG with PS, turn to the renormgroup FE g(xt, g) = g(x, g(t, g)). Its general 
solution is known; it depends on an arbitrary function of one argument - see Eq . (|19"D . 
However, at the moment, we are interested in a special solution linear in the second 
argument: g(x,g) = gX(x). The function X(x) should satisfy the equation X(xt) = 
X(x)X(t) with the solution X(x) = x v . Hence, g(x,t) = gx u . This means that in our 
special case, linear in g, the RG transformation (|T0| ) is reduced to the PS transformation, 

Rt =► S t : {x' = xt~\ g' = gt v } . (30) 

Generally, in RG, instead of a power law, we have an arbitrary functional dependence. 
Thus, one can consider transformations (0), ( |13"1) and ( |l6l) as functional generalizations 
of usual (i.e., power) self-similarity transformations. Hence, it is natural to refer to them 
as to the transformations of functional scaling or functional (self) similarity (FS) rather 
than to RG-transformations. In short, 

RG = FS , 

with FS standing for Functional Similarity^. 

Now we can answer the question on the physical meaning of the symmetry underlying 
FS and the Bogoliubov renormgroup. As we have mentioned, it is not a symmetry of a 
physical system or of equation(s) of the problem at hand, but a symmetry of a solution 
considered as a function of the relevant physical variables and suitable boundary param- 
eters. A symmetry like that can be related, in particular, to the invariance of a physical 
quantity described by this solution with respect to the way in which the boundary con- 
ditions are imposed. The changing of this way constitutes a group operation in the sense 
that the group composition law is related to the transitivity property of such changes. 

Homogeneity is an important feature of a physical system under consideration. How- 
ever, homogeneity can be violated in a discrete manner. Imagine that such a discrete 
violation is connected with a certain value of x , say, x = y. Here, RG transformation 
with the canonical parameter t has the form (|13|). 

The symmetry connected with FS is a very simple and frequently encountered property 
of physical solutions. It can easily be "discovered" in numerous problems of theoretical 
physics like classical mechanics, transfer theory, classical hydrodynamics, and so on 0, 
!3T| - see, above, Section 2.1.2. 



6 This notion was first mentioned in B9] and formally introduced [p0[ in the beginning of 80s. 
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3 Symmetry of solution in Mathematical Physics 

3.1 Constructing RG-symmetries and their use 

From the discussion in Sections 1.1 and 1.2 it follows that FS transformation in QFT is 
the scaling transformation of an independent variable x (and, possibly, the parameter y) 
accompanied by a functional transformation of the solution characteristic g. It is intro- 
duced by means of either finite transformations fllOD, (|T3|) and fllTf ) or the infinitesimal 
operator (0). Hence, the symmetry of a solution, i.e., FS symmetry, is commonly under- 
stood in QFT as the Lie point symmetry of a one-parameter transformation group defined 
by the operator of the (0)-type. 

Now, we are interested in getting answers to the following questions: 

• is it possible to extend the notion of RG symmetry (RGS) and generalize the form 
of RGS implementation that may differ from that given by (|7|) ? and if "yes" , 

• is it possible to create a regular algorithm of finding these symmetries ? 

The answer is positive to both the questions, and below we demonstrate the regular 
algorithm of constructing RGS in mathematical physics that up to now has been devised 
only for boundary value problem (BVP) for the (system of) differential equation(s) which 
we shall refer to as basic equations (BEs). The point is that these models can be ana- 
lyzed by methods of Lie group analysis which employ infinitesimal group transformations 
instead of the finite one. 

The general idea of the algorithm is to find a specific renormgroup manifold 1ZM. 
that contains the desired solution of BVP. Then construction of a RGS, that leaves this 
solution unaltered, is performed by standard methods of a group analysis of DEs. 

The regular algorithm of constructing RGS (and their application) can be formulated 
in a form of a scheme^] which comprises a few steps. It is illustrated in Figure 1. 

I. First of all, a specific renormgroup manifold 1ZM for the given BVP should be 
constructed which is identified below with a system of the fcth-order DEs 

F a (z,u,u {1) ,...,U( k )) = 0, a = l,...,s. (31) 

In ( pip and what follows we use terminology of group analysis and the notation of differ- 
ential algebra. In contrast with the mathematical analysis, where we usually deal with 
functions u a , a = 1, . . . ,m of independent variables x l , % — 1, . . . ,n and derivatives 
uf(x) = du a /dx\ ufj(x) = d 2 u a /dx l dxi , . . . that are also considered as functions of x, 
in differential algebra we treat u a , u", uf,, ... as variables as well. Therefore, in differen- 
tial algebra we deal with an infinite number of variables 

x = { x 1 }, u = {u a }, « (1 ) = {<}, M (2 ) = « i2 }, ••• ; (*,ii,... = l,...,n), (32) 



7 In the present form this scheme was described in pT|. One can find there historical comments and 
references on the pioneering publications. 
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Figure 1: RGS constructing and application to BVP in Mathematical Physics. 



where x % are called independent variables, u a dependent variables and um, Um, . . . deriva- 



tives. A locally analytic function f(x,u,u 



(i)' 



u 



(k), 



of variables 



with the highest 



fcth-order derivative involved, is called a differential function of order k. The set of all 
differential functions of a given order form a space of differential functions A, the universal 
space of modern group analysis [p^] -[p5|1. 

A particular way of realization of the first step can hardly be described uniquely, as 
it depends on both a form of basic equations and a boundary condition; generally, 1ZM. 
does not coincide with BEs. We indicate here a few feasible routines for this step. 

• One can use an extension of the space of variables involved in group transformations. 
These variables, e.g., may be parameters, p = {p^},j = l,...,l entering into a 
solution via the equations and/or boundary conditions. Adding parameters p to 
the list of independent variables z = {x, p} we treat BEs in this extended space as 
TZM. (|31]). Similarly, one can extend the space of differential variables by treating 
derivatives with respect to p as additional differential variables. 

• Another possibility employs reformulating of boundary conditions in terms of em- 
bedding equations or differential constraints which are then combined with BEs. The 
key idea here is to treat the solution of BVP as an analytic function of independent 
variables and boundary parameters b = {x 1 q ,Uq} as well. Differentiation with re- 
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spect to these parameters gives additional DEs (embedding equations) that, together 
with BEs, form 1ZM.. In some cases, while calculating Lie point RGS, the role of 
embedding equations can be played by differential constraints (for details see []3"T|) 
that come from an invariance condition for BEs with respect to the Lie-Bdcklund f\ 
symmetry group. 

• In the case when BEs contain a small parameter a, the desired 7£.A/f can be obtained 
by simplification of these equations and use of "perturbation methods of group 
analysis" (see Vol.3, Chapter 2, p. 31 in [J34fl). The main idea here is to consider 
a simplified (a = 0) model, which admits a wider symmetry group (see examples 
in the Section |4.2j below) in comparison with the case a^fl. When we take the 
contributions from small a into account, this symmetry is inherited by BEs, which 
results in the additional terms, corrections in powers of a, in the RGS generator. 

II. The next step consists in calculating the most general symmetry group Q that 
leaves the manifold 1ZM. unaltered. The term "symmetry group" , as used in the classical 
group analysis, means the property of the system (|3~T|) to admit a local Lie group of point 
transformations in the space A. 

The Lie algorithm of finding such symmetries consists in constructing tangent vector 
fields defined by the operator 

X = ed x < + T] a d u «, e,V a £A, (33) 

with the coordinates, £ 4 , n a that are functions of group variables and have to be deter- 
mined by a system of equations 

XF a =0, a=l,...,s, (34) 

(31) 

that follow from the invariance of 1ZM.. Here X is extended^ to all derivatives involved 
in F a and the symbol IQ means calculated on the frame fl3~T|). A system of linear ho- 
mogeneous PDEs (0) for coordinates £* , rj a , known as determining equations, is an 
overdetermined system as a rule. The solution of Eqs.fl3~3|) define a set of infinitesimal 
operators (^) (also known as group generators), which correspond to the admitted vector 
field and form a Lie algebra. In the case that the general element of this algebra 

X = E dXj , (35) 

3 

where A 1 are arbitrary constants, contains finite number of operators, 1 < j < I, the 
group is called finite- dimensional (or simply finite) with the dimension I; otherwise, for 

8 We use here the terminology adopted in Russian literature |52[ [35| . This symmetry is also known as 
generalized or higher-order symmetry |33| p4| . 

9 The extending of generators to the derivatives employs the prolongation formulas and is a regular 
procedure in group analysis (see, e.g. [p4| ) . 
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unlimited j or in the case that coordinates £\ 7] a depend upon arbitrary functions of 
group variables, the group is called infinite. 

The use of the infinitesimal criterion ( [51D for calculating the symmetry groups makes 
the whole procedure algorithmic and can be carried out not only "by hand" but using the 
symbolic packages of the computer algebra (see, e.g., Vol.3 in [J34J1 ) as well. In modern 



group analysis, different modifications of the classical Lie scheme are in use (see |32|, [33], [34 
and references therein). 

The generator (|33| ) of the group Q is equivalent to the canonical Lie-Backlund operator 

Y = K a d ua , K a = V a - f uf , (36) 

that is known as a canonical representation of X and plays an essential role in RGS 
constructing. 

However, the group defined by the generators (^) and ( |3"E| ) cannot yet be referred to 
as a renormgroup, as it is not related to a partial BVP solution of interest. 

III. To obtain RGS, the restriction of the group Q on particular BVP solution should 
be made which forms the third step. Mathematically, this procedure appears as checking 
the vanishing condition for the linear combination of coordinates k° of the canonical 
operator equivalent to (|3^) on a particular approximate (or exact) BVP solution U a (z) 






= . (37) 

u a = U a (z) 



Evaluating (37]) on a particular BVP solution U a (z) transforms the system of DEs for 
group invariants into algebraic relations^ Firstly, it gives relations between A 7 thus "com- 
bining" different coordinates of group generators Xj admitted by the 1ZM. ([!!]) . Secondly, 
it eliminates (partially or entirely) the arbitrariness that may appear in coordinates £ l , r] a 
in the case of an infinite group Q. In terms of the "classic" QFT RG terminology, where 
it exists only one operator X of ©-type (i.e., all Aj except one are equal to zero), the 
procedure of group restriction on a particular BVP solution g appr eliminates arbitrariness 
in the form of (3(g) -function. 

While the general form of the condition given by Eq.([J7]) is the same for any BVP 
solution, the way of realization of the restriction procedure in every particular case employs 
a particular perturbation approximation (PA) for the concrete BVP. 

Generally, the restriction procedure reduces the dimension of Q. It also "fits" boundary 
conditions into the operator (p5|) by a special choice of coefficients Aj and/or by choosing 
the particular form of arbitrary functions in coordinates £*, i] a . Hence, the general element 
(J35|) of the group Q after the fulfillment of a restriction procedure is expressed as a linear 

10 Similar relations were discussed in |32|| , Chapter 8, when constructing invariant solutions for the 
Cauchy problem for a quasi-linear system of first order PDEs. 
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combination of new generators Ri with the coordinates £*, r] a , 

X ^ R = J2 B J Rj , Rj = $8* + rjfd ua , (38) 

i 

where B-i are arbitrary constants. 

The set of RGS generators Ri each containing the desired BVP solution in its invari- 
ant manifold, define a group of transformations that we also refer to as renormgroup. 
Therefore, here we extend the notion of renormgroup and RG symmetry and the direct 
analogy with the "Bogoliubov RG" is preserved only for one-parameter group of point 
transformations. 

IV. The prescribed three steps entirely define the regular algorithm of RGS con- 
struction but do not touch on how a BVP solution is found. Hence, one more important, 
the fourth, step should be added. It consists in using RGS generators to find analytical 
expressions for the new, "improved" , solution of the BVP. 

Mathematically, this step makes use of RG=FS invariance conditions that are given 



by a combined system of fl31~|) and the vanishing condition for the linear combination of 
coordinates k" of the canonical operator equivalent to 



£#«? = £ £%-§<) =0. (39) 



One can see that conditions (|39|) are akin to (]37|) . However, in contrast with the previ- 
ous step, the differential variables u in ( |3"9"D should not be replaced by an approximate 
expression for the BVP solution U(z), but should be treated as usual dependent variables. 

For the one-parameter Lie point renormgroup, RG invariance conditions lead to the 
first order PDE that gives rise to the so-called group invariants (like invariant couplings 
in QFT) which arise as solutions of associated characteristic equations. 

A general solution of the BVP is now expressed in terms of these invariants. On the 
one hand, this is in direct analogy with the structure of RG invariant solutions in QFT - 
compare with Eqs. (p2|) and ([24|). 



On the other hand, it reminds the so-called n-theorem from the theory of dimensional 
analysis and similitude (see, Section 19 in [B21, Section 6 of Chapter 1 in [RB and historical 



comment to Section 43 in |3~7f ) directly related to power self-similarity, discussed above 



in Section 2.3. 



However, as we shall see later, in the general case of arbitrary RGS the group invariance 
condition obtained for BVP is not necessarily characteristic equations for the Lie point 
group operator. They may appear in a more complicated form, e.g., as a combination of 
PDE and higher order ODE (see Section |4l|). Nevertheless, the general idea of finding 
solution of the BVP as RG invariant solutions remains valid. 
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3.2 Examples of solution improving 

We present now a few examples of the RGS construction with further use of the symmetry 
for "improving" an approximate solution. 

3.2.1 Modified Burgers equation 

As the first example, we take the initial value problem for the modified Burgers equation 

u t - au 2 x - vu xx = , u(0, x) = f(x) . (40) 

It is connected with the heat equation 

Ut = VUxx (41) 

by transformation u = exp(aw/z/) and has an exact solution. Due to this, while using RGS 
to find a solution, one can check the validity of our approach. The RGS constructing for 
(|40|) is an apt illustration of the general scheme, shown in Figure 1 which may be helpful 
in understanding other examples of the general algorithm implementation. We review 
here in short the procedure and results of paper Ref. |&| . 



The RG-manifold 1ZM. (step I) is given by Eq.fl^OP with the parameters of nonlinearity 
a and dissipation v included in the list of independent variables. The Lie calculational 
algorithm applied to "KM. gives, for the admitted group Q (step II), nine independent 
terms in the general expression for the group generator 



X = ^ A 1 (a, v)X { + a{t, x, a, v)e~ aulv d u , (42) 

8=1 

Xx = Aut 2 d t + iutxd x - (u/a)(x 2 + 2vi)d u , X 2 = 2td t + xd x , 

X 3 = (l/u)dt , X A = 2vtd x - (v/a)xd u , X B = d x , X 6 = -{u/a)d u , 

X 7 = ad a + [(v/a) -u]d u , X s = 2vd v +xd x +2[u- {u/a)\ d u . 
Here, A 1 (a, v) are arbitrary functions of their arguments and a(t, x, a, v) is an arbitrary 



function of four variables, satisfying the heat equation (f4T|) . 

A set of operators A, form an eight-dimensional Lie algebra L§. The first six generators 
relate to the well-known symmetries of the modified (potential) Burgers equation (see, 



e.g., Vol.1, p. 183 in p4|). They describe projective transformation in the (t,x)— plane 
(Xi), dilatations in the same plane (A 2 ), translations along the t—, x— and u— axes (A 3 , 
A 5 and A 6 ) and Galilean transformations (A 4 ). The last two generators A7 and A 8 relate 
to dilatations of parameters a and v now involved in group transformations. 
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The procedure of restriction (step III) of the group (|42| ) admitted by 7\LA/f (|4T)| ) implies 
the check of the invariance condition fl37|) on a particular BVP solution u = U(t, x, a, v) 



\-q 00 + ^2A\a,v)K i \ = , Ki=r)i- £]u t - £fu x - t$u a - i\u v . (43) 

^ * =1 J u — U(t, x, a, v) 

This formula expresses the coordinate a of the last term in ( |4"2"D via the remaining coordi- 
nates of eight generators X; for arbitrary t, and hence for t = 0, when U(0, x, a, v) = f(x). 
As a result, we obtain the "initial" value a(0, x, a, v) and then, using the standard repre- 
sentation for the solution to the linear parabolic equation (f4"T|), the value of a at arbitrary 

8 

a(t, x, a, v) = — J^ A\a, v) < k^x, a, v) > . (44) 

Here, Kj(x, a, v) denote "partial" canonical coordinates Kj taken at t = and w = f{x). 
Symbol < F > designates the convolution of a function F with the fundamental solution 
of (f4"Tp, multiplied by the exponential function of / entering into the boundary condition 



< F(x, t, a, v) >= —rj= J dyF(y,t,a,v)exp 



{x-yf af(yY 



y/Airvt J \ Aid ' v I 

Substitution (|4]) in the general expression ( f42"l) gives the desired RG generators 

Ri = X % + Ql e- au l»d u , 

V 2 r 1 J-2 , n V 

Qi = - < x > , q 2 =< xf x > , g 3 = - < af x + vf xx > , q a = - < x > , 
a v a 

fe =< /x > , g 6 = - < 1 > > Qi =< J > , 6s = < xfx ~ 2/ + 2 - > . 

a a a 

Operators Ri form an eight- dimensional RG algebra RL 8 that has the same tensor of 
structural constants as L 8 , i.e. RL 8 and L s are isomorphic. Hence, the procedure of the 
group restriction eliminates the arbitrariness presented by the function a and "fits" the 
boundary conditions into RG generators by means of ft. 

It can be verified that the exact solution of the initial- value problem fl4"0|) 

is the invariant manifold for any of the above RGS operators. And vice versa, (|45|) 
can be reconstructed from an approximate solution with the help of any of the RGS 
operators or their linear combination. For example, two such operators, uR 3 = R t and 



[l/a)(R e + R 7 ) = R a were used in |38| to reconstruct the exact solution from perturbative 
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(in time and in nonlinearity parameter a) solutions. Below, we describe this procedure 
(step IV) using the operator R a , 

R a = d a + (l/o) (-u + e~ au ' v < f(x) >) d u . (46) 



It is evident that t, x and v are invariants of group transformations with ([46]), whilst finite 
RG transformations of the two remaining variables, a and u, are obtained by solving the 
Lie equations for (0), with £ the group parameter 

— = 1, a'| fc0 = a; — = a(t,x, a', i/) e aM '" - — , u'\ e=0 = u. (47) 

Combining these equations yields one more invariant J = e au / u — < 1 > for the RGS 
generator (|4*6]). Solution of (|4"T| ) along with ( PD gives the formulae for finite RG transfor- 
mations of the group variables {t, x, a, u, u] 

t' = t, x' = x , i/ = v, a' = a + £, u' = -^— In (e au/u + < e ef(x)/u - 1 >) . (48) 

a + £ v 7 

Choosing the value a equal to zero, which is a starting point of PA in a, we get a' = £. 
Then after excluding t, x, v and £ from the expression for v! ( |Hj|) and omitting accents 
over £', x 1 ,1/ , u' and a' the desired BVP solution (fiSD is obtained. It also follows directly 
from J7 in view of the initial condition J | a= o= 0. 

A similar procedure can be fulfilled for the other RG operator, 

Rt = d t + e~ aulv < af 2 x + uf xx > 8 U , 

which is consistent with the PA in time t. Although invariants for R t and finite RG 
transformations differ from that for fl46|), the final result, i.e., the exact solution of BVP 
(fin) given by (ff§), is the same. This possibility is the distinct demonstration of the 
multi-dimensional RGS to reconstruct the unique BVP solution from different PA: either 
in parameter a or in t (though we used only two one-dimensional subalgebras herd3). 

3.2.2 BVP for ODE: simple example 

Quite recently, the QFT renormalization group ideology has been applied, a bit straight- 



forward, in mathematical physics for asymptotic analysis of solutions to DEs [[11], |2| and 
in constructing an envelope of the family of solutions ||43|| . 

Our second methodological example with linear ODE is presented here in order to 
illustrate the difference between our approach and the "perturbative RG theory" devised 
in 



11 1 for a global analysis of BVP solutions in mathematical physics. 



Consider a linear second order ODE for y(t) with the initial conditions at t = r, 

y tt + yt + ey = 0, v(t)=u, y t (r)=w, (49) 



11 This can be considered as a construction parallel to the one used in Ref. [p9[ 
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which has the exact solution: 

V = C + e-»<t-r) + CLe-r-e-) , 7± = ^ , K = VT^Te , C ± = T ^±^ . (50 ) 

Provided that the parameter e is small, the solutions to Eq.(f49|) has been treated 
in [41 1 with the aim to demonstrate effectiveness of the "perturbative RG theory" for 
an asymptotic analysis of a solution behaviour. The main goal of this treatment was 
to improve a perturbative expansion in powers of e with secular terms oc e(t — r) and 
obtain^ a uniformly valid asymptotic of a solution 

y = c +e (- 1+£ ( 1+£ » (^ + c_e- £ ( 1+£ ) ^ + 0(e 2 ) , (51) 

c+ « - ((1 + 2e)w + eu) , c_ « ((1 + 2 £ )w + (1 + e)fi) , 

which is accurate for small values e <C 1 but for arbitrary values of the product e(t — r). 

We are going to show that the use of our regular RG algorithm enables one to improve 
a PA solution (either in powers of e or in t — r) up to the exact BVP solution (pTf). 



Rewriting fl49|) in the form of the system of two first order ODEs for functions u = y 
and w = y t , 

u t = w , Wt — —ew — u , (52) 

we construct 1ZM (step I) using the invariant embedding method (this approach has 
first been realized in [|44| ). Then, 7£.M is presented as a joint system of BEs (E^) and 
embedding equations 

m t — (eiy + -u)m^ — wUu = , w r — (eu) + u)ww — wwn = , 

treated in the extended space of group variables which include the parameters r, w, u of 
boundary conditions in addition to t and dependent variables u, w. 

Omitting tedious calculations related to the following two steps (steps II and III ), 
we present here two examples of resulting RGS generators 

R T = d T — (w + eu)du, + wda , 

Re = d £ - (\{2w + u)(l - t/2) + (t/2)uj w + \ {2w + u) d u 

- (\{2w + u)(l - t/2) + (r/2)uj ^ + -^ (2t5 + u) 9 S , 



12 The algorithm used in |4l| for improving PA solutions with secular terms involves a) an introduction 
of some additional parameters in solutions, b) a special choice of these parameters that eliminates secular 
divergencies, and c) imposing a condition of independence of a solution upon the way of introducing 
these parameters. In some cases, this algorithm, directly borrowed from QFT RG-method, gives an exact 
solution. However, the question of correspondence of this construction to a transformation group of a 
solution of BEs remains open. 
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that involve the initial values w, u and initial point r in RG transformations. In addition, 
R e transforms the parameter e. 

Now, the procedure of constructing the BVP solution (EDf) (step IV) is similar to that 



used in the previous Section |3.2.1| and employ finite transformations that are defined by 
the Lie equations for the operators (|K|). For R T functions u, w and the parameter e 
are group invariants, while the translations of r and the corresponding transformations 



of u, w restores the exact solution (50) from the PA in powers of t — r (note that the 
parameter e is not necessarily small in this PA!). 

For R £ the difference t — r is group invariant, whilst the transformation of e and related 
transformations of u, w, u, w restore the exact solution ( |5UD from the PA (discussed in |¥T| ) 
in powers of e. Hence, as in the previous Section [3.2.i| , both the RGS generators (|53|) 
reconstruct the unique BVP solution from different PAs. 

4 RG in Nonlinear optics 

4.1 Formulation of a problem 

As a problem of real physical interest, take BVP that describes self- focusing of a high- 
power light beam. While the problem plays an important role in nonlinear electrodynam- 
ics since 60s, the detailed quantitative understanding of self-focusing is still missing |45| , 



and there is no method which allows to find an analytic solution to the corresponding 
equations with arbitrary boundary conditions. 

Here, we demonstrate the great potential of the RGS approach in constructing ana- 
lytic solutions of BVP equations with arbitrary boundary conditions. The RGS method 
allows to consider different types of BEs for self-focusing process which include plane and 
cylindrical beam geometry, nonlinear refraction and diffraction. The merit of the RGS 
method is that it describes BVP solutions with one- or two-dimensional singularities in 
the entire range of variables from the boundary up to the singularity point. 

Let us start with BVP for the system of two DEs 

v z + vv x — an x = , n z + nv x + vn x + [y — 1) {nv/x) = , (54) 

v(0,x) = 0, n(0,x) = N(x), (55) 

which are used in nonlinear optics of self-focusing wave beams when diffraction is negli- 
gible. 

We study spatial evolution of the derivative of the beam eikonal v and the beam 
intensity n in the direction inwards the medium z and in the transverse direction x. 
The term proportional to a is related to nonlinear refraction effects; v = 1 and v = 2 
refers to the plane and cylindrical beam geometry, respectively. Boundary conditions 
(|55|) correspond to the plane front of the beam and the arbitrary transverse intensity 
distribution. 
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4.2 Plane geometry 

In the plane beam geometry (at v = 1) Eqs. (|54j) can be reduced to the system of BEs 

r w - nx„ = , Xw + ar n = , (56) 

for functions t = nz and x = z — vz of w = f/a and n arguments, with boundary 
conditions 

r(0,n)=0, x(0,n) = fT(n), (57) 

where H(n) is the inverse to N(x). Here, the procedure of RGS constructing makes 



use of the Lie-Backlund symmetry and is described as follows |H]]. The manifold 7<LA/f 



(step I) is defined by Eqs. (|56|) treated in the extended space that include dependent and 
independent variables r, x, w, n and derivatives of r and x with respect to n of an 
arbitrary high order. The admitted symmetry group Q (step II) is represented by the 
canonical Lie-Backlund operator 

X = fd T + gd x , (58) 

with the coordinates / and g that are linear combinations of r and x and their derivatives 
d % rjdn % and d l x/dn l , % > 1 with the coefficients depending on u> and n. 

The restriction of the group admitted by 1ZM. ( |56|) (step III)) implies the check of 
the invariance condition (|37|) that yields two relations 



/ = 0, y = 0. (59) 

These relations should be valid on a particular solution of BVP with the boundary data 
([57]). For example, choosing the so-called "soliton" profile, N(x) = cosh~ 2 (x), i.e., H{n) = 



Arccosh (l/y/n), we have 

/ = 2n(l - n)r n „ - nr n - 2nw(x n + nx nn ) + (aw 2 / 2) nr nn , 

(60) 
g = 2n(l - n)xn« + (2 - 3n)x„ + aw (2nr nn + r n ) + (aw 2 / 2) (nx« n + Xn) ■ 

Dependence on r nn and Xnn indicates that here RGS is the second-order Lie-Backlund 
symmetry. In order to find a particular solution of a BVP (step IV), one should solve 
the joint system of BEs fl5"6|) and second-order ODEs that follow from the RG=FS invari- 



ance conditions (|59D and (|60|) . The resulting expressions [46|] - the well-known Khokhlov 
solutions^ 

v = —2anz tanh(x — vz) , an 2 z 2 = ncosh (x — vz) — 1 , (61) 

describe the process of self-focusing of a soliton beam: the sharpening of the beam inten- 
sity profile with the increase of z is accompanied by the intensity growth on the beam 



13 



In Ref . Jl7| , where this solution was first obtained, it did not result from a regular procedure. 
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axis. The solution fl6"T|) is valid up to the singularity point where the derivatives v x and 
n x tend to infinity whilst the beam intensity n remains finite 

<4 = l/2^, n s ° l ng = 2. (62) 

Here, the Lie-Backlund RGS enables one to reconstruct the BVP solution and describe 
the solution singularity for the light beam with the soliton initial intensity profile. One 
more example of an exact BVP solution obtained with the help of Lie-Backlund RGS 
(with the initial beam profile in the form of a " smoothed" step) can be found in |36|] . 

For arbitrary boundary data, it turns to be impossible to fulfill the condition (|59|) 
with the help of the Lie-Backlund symmetries of any finite order, and one is forced to 



use a different algorithm |31|, |46| of RGS constructing, based on the approximate group 
methods. Here, (step I)) 1ZM. is given by BEs ( |56|) with a small parameter a, and 
coordinates of the group generator (^) (and, hence, coordinates of the RGS operators) 
appear as infinite series in powers of a 

oo oo 

/ = X>T; = 5>y. ( 63 ) 

j=0 4=0 

The procedure of finding the coefficients f l , g l (step II)) leads to the system of recurrent 
relations that express higher-order coefficients f l+1 , g l+l in terms of previous ones f l , g l . 
It means that once the zero-order terms are specified, the other terms are reconstructed 
by the recurrent relations. 

The coefficients f l and g % contain an arbitrary function of n, x\s] an d Tr s i — w(sx[ s ] + 
n%r s+ i]) where subscript [s] denotes the partial derivative of the order s with respect to 
n. This arbitrariness is eliminated by the procedure of group restriction (step III)), i.e., 
by imposing the invariance condition fl59[). For particular forms of /° and g°, that is 
for partial boundary conditions ([57]), infinite series are truncated automatically, and we 
arrive at the exact RGS. 

One example of this kind is given by Eqs.([K]) that have a binomial structure / = 
f° + af 1 , g = g° + ag 1 . If we neglect the higher-order terms in the case of arbitrary 
boundary conditions (when series (^) are not truncated automatically), then we get an 
approximate RGS which produces an approximate solution to the BVP. As an example, 
we give here two sets of expressions for the coordinates f l and g l for the Gaussian initial 
profile with N(x) = exp(— x 2 ), i.e., H{n) = (ln(l/r7,)) ' , which define approximate RGS 

a) /° = 1 + 2n X Xn , 9° = , f 1 = -2rr n + — , g 1 = -2 {r X n + xr n ) , 

n 

b) f° = 2n(rxn + r n x) , g° = 1 + 2n\Xn , f 1 = ^X^a , Q 1 = 2 (xXa ~ rr n ) . 

Here, linear dependencies of / and g upon first-order derivatives indicate that RGS 
is equivalent to Lie point symmetry. The peculiarity of the case b is a dependence of 
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/ and g not only on derivatives with respect to n but also with respect to a: it means 
that the parameter a is also involved in group transformations. In the non-canonical 
representation (PBf), the RGS generator in this case has the form 



Rcuassi = 2rd w + 2n\d n + 2a\d a - d x . (64) 

The last step IV is performed in a usual way by solving the joint system of BEs 
(|56|) and equations that follow from the RG=FS invariance condition ( (j9|) , or else, using 
invariants of associated characteristic equations for RG operator provided that RGS is a 
Lie point symmetry. We give here the solution that follows from RGS (164]), 



2 (■, r, 2\ 2 1 1 2xanz ,„_. 

x = (1 — lanz hi—; — , v— -. (65) 

v ' nil — anz 1 ) 1 — 2anz 2 

These expressions describe a self-focusing Gaussian beam (the plot n(x) for this solution 
is presented at the end of the section on Figure 2), that is qualitatively very similar to 
the spatial evolution of the soliton beam fl6l|) . Moreover, the singularity position and the 
value of maximum beam intensity at this point coincide with analogous values ( |6"2"D for the 
soliton beam. Although formulae fl65| ) correspond to an approximate BVP solution, they 
exactly describe the behaviour of n on the beam axis at x — 0. To estimate the reliability of 
result fl65|) in the off-axis region, we compared it with another approximate BVP solution 
which arises from the approximate RGS in the case a. These approximations agree very 
well (details are presented in fJ6||), thus proving the accuracy^ of the RG approach. 

4.3 Cylindrical geometry 

In the above discussion we dealt with the plane beam geometry and took into account 
only effects of nonlinear beam refraction, neglecting diffraction. The flexibility of RGS 
algorithm allows one to apply it in a similar way to a more complicated model as compared 



to (|56"D, e.g., for the cylindrical beam geometry, v — 2. Omitting technical details, we 
present the RGS generator for the cylindrical parabolic beam with N = 1 — x 2 



Rpar = { 1 — 2az ) d z — 2azxd x — 2a(x — vz) d v + Aanzd n 



The BVP solution is expressed in terms of group invariants for this generator: 

2 
■' - ~ ~ 2 2 , XV f-i <-, 2 



J\ = — ; 3i = ng; J 3 = 2ax - v g + —g z ; g = (1 - 2az j . (66) 

g l 

The explicit form of dependencies of J% — 1 — Jf , Cfz = 2a Ji upon J7i follows from 



the boundary conditions (|5q). They lead to the well-known solution [47 



v = {x/{2g))g z , n = {l/g){\ - {x 2 / g)) (67) 



14 One more evidence is provided by the comparison of approximate and exact BVP solution for the 
soliton beam performed in Bq]. 
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that describes the convergence of the beam to the singularity point z p s ^ g = l/v2a where 
g — and n —>■ oo. The solution singularity is two-dimensional here: the infinite growth 
of beam intensity in the vicinity of the singularity z — > z v °^ g is accompanied by the infinite 
growth of the derivative v x and collapsing of the beam size in the transverse direction. 

The RGS algorithm based on approximate group methods can also be applied in 
the case when besides nonlinear refraction also diffraction effects are taken into account. 



Then, the first equation in (|54 ) should be modified by adding the diffraction term 



-(3d x {(x l -»/^) d x (aT-^V^)} 



Standard calculations done in compliance with a general scheme for thus modified 1ZM. 
(for details see ||48| ) give the RGS generator for the cylindrical beam geometry [y = 2) 



" : ' 1 + "7; Sxx + ^ Sx 



d n - 

(68) 



RGauss2 = (l + Z 2 S X x) d * + ( zS X + V^^XXj 9 X + S X^v ~ 

Here the function S, defined by the form of the intensity boundary distribution, 

S(x) = aN(x) + — 7=T d * (x d xjNtx)) , 

contains two small parameters, a and (3, and, as in the case (3 = 0, there exist specific 
forms of boundary distribution, N, for which the RGS operator (p8|) defines exact (not 
approximate) symmetry valid for arbitrary values of a and /3. 

Constructing a particular BVP solution (step IV) implies the use of group invariants 
related to (p8|), and the procedure is similar to that one for the parabolic beam. For the 
Gaussian wave beam, N = exp(— x 2 ), the result is as follows: 

x — x 2 ^(/3- ae~ x J 
v(z,x) = , n(z, x) = e~ M -- f- -^- . (69) 

Here \ and \i are expressed in terms of t and x by the implicit relations 

/3 (/i 2 - x 2 ) + a (e^ 2 - e^ 2 ) = 2^V (/3 - 

X (l + 2z 2 (/3 






x = y 1 -+• Z2: d — ae 



The solution ( |5"9"D describes the self-focusing of the cylindrical Gaussian beam that 
gives rise to the two-dimensional singularity: both the beam intensity n and derivatives 
v x ,n x go to infinity at the point z^°^ g ss = 1/ J2{a — (3) provided that a > (3. A detailed 
analysis of (|5"9"D and more general solutions with a parabolic form of an eikonal at z = 0, 



v(0, x) = —x/T, is given in |48|, E^j 
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Figure 2: Intensity n versus transverse coordinate x for a plane (left panel) and cylindrical 
(right panel) beam geometry for a few values of distance z from the boundary z 2 > Z\ > 
z = 0. 



To illustrate the difference between the one- and two-dimensional solution singular- 
ities, in Figure [| we present a typical behavior of the wave beam intensity, defined by 
Eq.(|B5D and floTf). The left panel corresponds to the plane beam geometry, v — 1, and 
without diffraction, /3 = 0, while the right one is concerned with a cylindrical wave beam, 
v = 2, with both nonlinearity and diffraction effects included. Diverse curves describe 
beam intensity distribution upon coordinate x at different distances from the medium 
boundary, where we have the collimated Gaussian beam, N = exp(— x 2 ). 

It is clear that in the plane geometry the derivative of the beam intensity with respect 
to x turns to infinity at some singular point, while the value of intensity on the axis 
remains finite. In cylindrical case the solution singularity is two-dimensional: both the 
beam intensity and its derivative with respect to transverse coordinate turn to infinity 
simultaneously at point z S i ng . This last example demonstrates the possibility of the RGS 
approach to analyse two-dimensional singularity. In the practice of RG application to 
critical phenomena, this correlates with the case of "two renormalization groups" ™' ! 



5 Overview 

To complete our review, we indicate milestones in evolution of the RG concept. Since its 
appearance in QFT, RG served as a powerful tool of analyzing diverse physical problems 
and improving solution singularities disturbed by perturbation approximation. 

The development of the RG concept can be divided into two stages. The first one (since 
the mid-50s up to the mid-80s) is summarized in Figure 3. Besides early history (discovery 
of RG, formulation of the RG method and application to UV and IR asymptotics), it 
comprises the devising of the Kadanoff- Wilson RG in the 70s and following explosive 
expansion into other fields of theoretical physics. 

During this stage, the formulation of the RG method was based on the unified scaling 
transformation of an independent variable (and/or some parameters) accompanied by a 
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Figure 3: Early development of concept: from Bogoliubov RG to Wilson RG and FS. 



more complicated transformation of a solution characteristic g^ = g(fi,g) - see Eqs.fllOD, 
(0) and (|g) in Section [O. 

Here, the main role of RG=FS was the a priori establishment of the fact that the 
solution under consideration admits functional transformations that form a group. 

Any particular implementation of the RG symmetry differs in the form of the func- 
tion^) g(n, g) (or (3(g)) which, in an every partial case, is obtained from some approximate 
solution. 

The next stage, after the mid-80s, is depicted in the Figure 4. The scheme describes 
the entire evolution of the Bogoliubov RG. There were several important reasons in further 
devising the RG concept in theoretical physics at this period. On the one hand, it was due 
to the extension of the notion of FS and RG symmetry that until then were based on one- 
parameter Lie group of point transformations. Appending multi-dimensional Lie point 
groups and Lie-Backlund groups to possible realization of group symmetry enhanced the 
capability of RG method. On the other hand, this additional possibility arose due to 
the mathematical apparatus that was used in mathematical physics to reveal RGS. The 
advantage came from infinitesimal transformations that enabled to describe RGS by an 
algebra of RG generators. However, in contrast to the situation typical of QFT models 
with only one operator, in mathematical physics we have finite or infinite-dimensional 
algebras. Both their dimension and the method of construction depend upon a model 
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Figure 4: Evolution of concept: from the Bogoliubov RG via FS to RG-symmetries. 

employed and upon a form of boundary conditions. 

The use of infinitesimal approach results in constructing the RG-type symmetry with 
the help of regular methods of group analysis of DEs. Precisely, this regular algorithm 
naturally includes the RG=FS invariance condition in the general scheme of constructing 
and application of RGS generators (see also our recent review ||9|]). Within the infinitesi- 



mal approach this condition is formulated in terms of vanishing of canonical RG operator 
coordinates, which is especially important for Lie-Backlund RGS because finite transfor- 
mations in this case are expressed as formal series. In particular, this property attribute 
a new feature to the RG analysis of a BVP solution with singular behavior, making a 
singularity analysis more powerful. 

At the same time, as the group analysis technique is still developing - here we mean 
both extension to new types of symmetries and application to more complicated mathe- 
matical models, e.g., including integro-differential equations - we have a clear perspective 
that the possibilities of a regular scheme based upon the Bogoliubov renormalization group 
method are far from being exhausted. 
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